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ABSTRACT 


Time-channel response functions have been derived for a time-of- 
flight spectrometer used to measure monoenergetic secondary neutrons in 
the energy region from 2 to 100 MeV. produced in homogeneous slab targets 
from the interactions of monoenergetic charged particles with the target 
nuclei. The effects on the flight-time distribution of the target and 
detector dimensions and orientation, the spatial distribution of the beam, 
and the pulse-height dependent time slewing (walk) were demonstrated. 
Response functions were derived assuming that the neutron emission from 
the target is independent of the azimuthal scattering angle and the polar 
angle over the angular region subtended by the detector, that neutron 
detection may be treated in the first collision approximation with the 
light being emitted uniformly from the interaction point and collected 
uniformly from all points in the detector volume, that the beam dimension 
perpendicular to the scattering plane can be ignored, and that the 
instrumental timing jitter can be approximated by a Gaussian distribution 
independent of the pulse height. The time-channel response functions 
were obtained from the integral over the time -analyzer response function 
of the convolution of the physical flight-time distribution of detected 
neutrons and the time-slewing response function. The maximum flight-time 
uncertainties, At/t, were derived in terms of the target and detector 
dimensions assuming that all dimensions have negligible effect except the 
one being evaluated. Comparisons are made between the computed distri- 
butions and those obtained analytically. The expected dominant effects 


xv 



of the detector thickness are demonstrated along with the dependence of 
the target dimensions on the neutron velocity and the angular difference 
between the neutron scattering angle and the target rotation angle . The 
distributions arising from the target and detector thickness were shown 
to be square except for attenuation effects, while those obtained for 
large detector radii showed a l/t dependence. The time distribution for 
large beam size reflects the beam spatial distribution. The numerical 
methods are described, and typical distributions illustrate each con- 
tribution to the overall response function. Where multiple scattering 
does not dominate, the time integral of the computed response functions 
are in good agreement with known detector efficiencies. 


CHAPTER I 


INTRODUCTION 

The determination of neutron energy by time -of -flight is one of 
the most direct methods for neutron spectroscopy. The method was origi- 
nally conceived for measurements of neutrons in the electron volt energy 
region, hut with the development of fast, highly efficient neutron 
detectors and improvements in pulse electronic instrumentation, measure- 
ments were extended into the region up to 20 MeV. 1-4 Recently, flight- 
time techniques have been used at even higher energies to measure target 
yields and differential cross sections for the interactions of 50- to 
150-MeV. protons with nuclei. 5 

Experiments have been performed at the Harvard Synchrocyclotron 
using time -of -flight techniques to measure the differential cross sec- 
tions for the production of neutrons resulting from the interactions of 
160 -MeV . protons with target nuclei. 6 As part of the analysis of these 
data, this study was initiated to determine the time response function 

1 F . G. J. Perey, Inelastic Scattering of 14-MeV. Neutrons in 

Carbon, Oxygen, and Lithium, Ph.D. dissertation. University of Montreal , 

I960. 

2 G. K. O'Neill, Phys. Rev . 1235 (1954). 

3 A. Adam, G. Palla, and P. Quittner, Acta. Phys. Hung. Tom. XVII 
Ease. 3> 1964. 

4 J. Rethneiv et al ., Nucl. Inst. Methods 17, 273 (1962). 

5 P. H. Bowen et al ., Nucl. Phys . 30, 475 (1962). 

6 R. W. Peelle et al.. Neutron Phys. Div. Ann. Progr. Rept. 

Sept. 1, 1962, 0RNL-3360, p. 286. 
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for monoenergetic neutrons in the range 2 to 100 MeV. in the experimental 
geometry . 

I. A SIMPLIFIED FLIGHT -TIME SPECTROMETER 

Figure 1 shows a simplified time -of -flight spectrometer based on 
the one described by Peelle et al . 7 The flight-time of the neutron is 
inferred from the time interval between the pulse produced in the zero- 
time detector by the incident proton and the pulse from the disc -shaped 
NE -21^ neutron detector. (Not shown in the diagram are the additional 
detectors and instrumentation for charged-particle and gamma-ray 
discrimination since it is assumed that only neutrons interact in the 
neutron counter.) A scintillating detector is shown as the zero-time 
detector; however, such devices as induction electrodes or the rf signal 
from the proton accelerator may be used. In some applications, the alpha 
or gamma radiation accompanying the formation of the neutron is detected 
by a counter located near the target to provide the zero-time signal. 

The conversion gain of the time analyzer is 0.500 nsec ./channel . 

II. REVIEW OF NEUTRON FLIGHT-TIME SPECTROSCOPY 

The resolution of a time-of-flight spectrometer is determined for 
a given neutron energy from the flight -time interval measurement 

7 R. W. Peelle et al .. Differential Cross Sections for the Produc- 
tion of Protons in the Reactions of 160-MeV Protons on Complex Nuclei , 
ORNL-3tt&7> September 1966. 

S NE-213 is a liquid scintillator manufactured by Nuclear 
Enterprises, Ltd., Winnipeg, Canada. 
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uncertainty combined with the uncertainty in the flight path . 9 The 
timing uncertainty is the combination of the partial uncertainties 
arising from the scintillator detector, photomultiplier tubes, timing 
circuits, and time analyzer channel width. The flight-path uncertainty 
is introduced by the variation in the birth and detection points within 
the target and detector. The relative uncertainty in the measurement of 
the energy is twice as large as the combination of all geometric and 
timing uncertainties. To achieve optimum energy resolution, the dominant 
contributions must be controlled. 

A detailed discussion of time-of -flight techniques is given by 
Neiler and Good . 10 They separate the instrumental resolving time into 
three "measurable" functions; the uncertainty in the time -of -origin of 
the neutron, the geometry of the detecting medium, and the uncertainties 
associated with detection and electronic processes. Each of these fac- 
tors requires special consideration by the experimenter. Rybakov and 
Siderov 11 point out the improved resolving power is, in general, obtained 
by increasing the flight path. This idealization is not always practical 
since the loss in counting rate with distance and the signal to back- 
ground ratio for available beam intensity may require short flight paths. 
Additionally, scattering from the floor and walls of the experimental area 

9 The uncertainty is defined here as the predictable spread in 
values for a given measurement and not the unknown systematic error in 
the mean of the resulting distribution. 

10 J. H. Neiler and W. M. Good, Fast Neutron Physics, (Interscience 
Publishers, New York, 1963 ), pp. 509-621. 

11 B. V. Rybakov and V. A. Siderov, Fast Neutron Spectroscopy, 
(Consultants Bureau, Inc., New York, i960), pp. 101 - 109 • 
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requires more complicated apparatus when using long flight paths, because 
of the possibility of detecting neutrons degraded in energy. 

Brooks 12 lists four criteria for selecting a detector for flight- 
time measurements. These are: time-resolution, detector efficiency, 

discrimination against background, and physical considerations, e. g., 
thickness and/or area. The timing resolution that can be achieved depends 
on the choice of the scintillator and photomultiplier tube. The large 
amplitude range in the light output from the detection event causes a 
portion of the amplitude spread to be translated into time spread by the 
electronics normally called time walk or slewing. Additionally, the 
statistical nature of the cathode photoelectron emission and electron 
multiplier introduce timing fluctuation or jitter. As Lundby 13 has 
shown, the construction of the scintillator-phototube assembly must also 
be considered. Detailed studies of the timing resolution have been made 
by Colombo, Gatti, and Pignanelli 14 who have related the resolving time 
to the variance of the centroid of the current pulse from the phototube 
and the pulse-height dependent time-slewing. Graphic results have also 
been presented by Gatti and Svelto 15 for estimating the time resolution 
for given bias levels as a function of the scintillator decay time 
constant and statistical properties of the phototube specified by the 

12 F. D. Brooks, Proc . Symp. on Neutron Time -of -Flight Methods, 
Saclay (1961), p. k03 . 

13 A. Lundby, Rev. Sci. Instr . 22, 32^. (1951). 

14 S. Colombo, F. Gatti, and M. Pignanelli, Nuovo Cimento 5, 1739 

(1957). 

ls E. Gatti and V. Svelto, Nucl. Inst, and Methods 30j 213 (196^-). 
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transit time fluctuations. Several methods of time measurement and 
their associated relative timing errors are discussed in detail by 
Neiler and Good and Rybakov and Siderov. 

At neutron energies over 0.5 MeV., detection is usually accom- 
plished with n-p scattering using organic scintillators. Organic scintil- 
lators are characterized by fast decay times and relatively high effi- 
ciency. Birks 16 has investigated the basic processes of scintillation 
and has derived analytic relationships between the fast and slow light 
components of the decay process while Batchelor et al . 17 and Verbinski 
et al. 1 8 have separately measured the response functions of organic 
scintillators for fast neutrons and have established the relationships 
between the light output and the energy for electrons, protons, alpha 
particles, and carbon recoils. Analytic methods for determining the de- 
tection efficiency of organic phosphors to neutrons in the energy region 
below 20 MeV. have been suggested by Kurz 19 and by Schuttler 20 and abso- 
lute efficiency values have been reported by Love and his coworkers. 21 

16 J. B. Birks, IRE Transactions on Nuclear Science, NS^£, 2 (i960). 

17 R. Batchelor et al ., Nucl. Instr. Methods 15, 20 (1962). 

18 V. V. Verbinski et al .. The Response of Some Organic Scintillators 
to Fast Neutrons, Shielding Division Report, ANS-SD-2, December 19&T. ~~ 

19 R. J. Kurz, A 709/7090 Fortran II Program to Compute the Neutron- 
Detection Efficiency of Plastic Scintillator for Neu tron Energies From 
1 tO 500 MeV ., UCRL-11539, March 1964. 

20 R. S. Schuttler, Efficiency of Organic Scintillators to East 
Neutrons, ORNL-3888, July 1966. 

21 T. A. Love et al .. Absolute Efficiency Measurements of NE-215 
Organic Phosphors for Detecting lk.k- and 2.6-MeV. Neutrons, 0RNL-3893, 
September 1966. 
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Recalling the detector criteria of Brooks, one must finally consider 
the physical dimensions of the detector. For thin detectors in which the 
attenuation is minimal, the efficiency is proportional to the detector 
thickness. For detectors of thicknesses comparable to the fast neutron 
mean free path, the proportionality breaks down, e. g., 7 cm. at 2 MeV. 
Grismore and Parkinson 22 have investigated the relationship between the 
resolving power and the neutron energy for flight paths' of 5-3 and 15 m. 
and detector thicknesses of 2 and 8 cm. They have shown that for 6-MeV. 
neutrons produced by 16 0(d,n) reactions and 3 nsec, time resolution, the 
geometric factors which limit the resolving power are just equal to the 
time resolution when a 6-cm. detector and a 5 -cm. thick gas target are 
used. Rybakov and Siderov give graphic results for the neutron energy 
and flight distance required to obtained 3 °jo and 1 $ energy resolution for 
a 5 -cm. thick detector for At = 3 nsec. 

Typical factors which determine the resolution time have been 
summarized in a paper by Smith 23 for a time -of -flight spectrometer used 
to study neutrons with energies between 2 and 4 MeV. 

III. PURPOSE OF THE THESIS AND PROPOSED WORK 

The purpose of this thesis is to calculate the expected neutron 
flight-time distributions for monoenergetic neutrons resulting from the 
bombardment of homogeneous targets by monoenergetic protons and to study 

22 R. Grismore and W. C. Parkinson, Rev. Sci. Instr . 28, 2^5 (1957)* 

23 R. V. Smith, Conference on Neutron Physics by Time -of -Flight , 
Gatlinburg, Tennessee, 1958, p. 103* ” 
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the behavior of these distributions for variations in the geometric 
parameters of the beam, target, and detector. 

Several restrictions have, however, been imposed on the design of 
the model spectrometer which is used as well as the physical phenomena 
leading to birth and detection of the neutrons : 

1. The target and detector are right -circular cylinders mounted so 

that the particles first strike the flat surfaces . 

2. The calculations are made for 160-MeV. protons whose specific 

energy loss in the target is approximated by its true value 
at the center of the target. Multiple scattering of the protons 
or secondaiy neutrons in the target is ignored. The proton 
density distribution is approximated either by a bivariate normal 
frequency function with a diameter much less than the target 
diameter, or as uniform over the target . 

3 . Neutron emission from the target is assumed independent of the 

asimuthal scattering angle, as well as independent of the polar 
angle over the angular range subtended by the detector. No 
interactions occur in the region between the target and detector, 
k. Neutron detection is restricted to the first collision approxima- 
tion and the resulting light is emitted from the interaction 
point and collected uniformly from all points in the detector 
volume. Light transit times are obtained from the normal dis- 
tance to a flat photocathode with the scintillator coupled 
directly to the multiplier tube. Phototube delays are taken as 
independent of the position where the light strikes the photo- 


cathode . 
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5- Timing jitter was represented by a normal frequency function with 
standard deviation independent of pulse height from the 
scintillation detector. 

For various secondary neutron energies the analysis will show the 
effects on the flight -time distribution of: 

1. the angle between the beam direction and the target -detector axis, 

2. the target thickness and orientation, including consideration of 

the energy loss of the incident protons, 

3. the detector dimensions, 

4. the spatial distribution of the proton beam, 

5. the timing fluctuations between the flight time and the observed 

"machine " time } and 

6. the time-slewing and analyzer response. 


CHAPTER II 


COORDINATE SYSTEMS 
I . SPACE COORDINATES 

Figure 2 shows the target and detector in the typical experimental 
configuration treated in the calculations which follow. A cartesian 
coordinate system, x-y-z, is constructed such that the positive direction 
of the x-axis is along the direction of the proton beam with the x-y plane 
containing the center of the detector and the target normal. The target 
coordinate system, £-r)-Q> is constructed so that the vertex of the £-tj-£ 
system is in conjunction with the origin of x-y-z with £ and 77 constrained 
to rotate in the x-y plane and £ along the target normal. The angle of 
rotation, Y, is limited to values defined by |y| < n/2. The target 
thickness, T, is measured along | such that r - = T and 

T 1 £ {S T 2 II-l 

for all values of 4 in the target . 

The position of the neutron detector is given by Id s L and 0, 
where L is a vector in the x-y plane from the origin to the face of the 
detector and normal to it, and Q is the angle between the x-axis and L. 

The detector coordinate system, x' -y' -z' , is constructed within 

— > 

the detector so that the origin of x' -y' -z' is at the terminus of L with 
the positive direction of the x'-axis along the axis of the detector. 
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Transformation from the primed system to the unprimed system is accomplished 
using the expressions 


x = x' cos 6 - y' sin0 + |l| 

y ~ y' cos 0 + x' sine II-2 

z = z' 

The distance from where the proton intersects the face of the 
target to the point of interaction in the target, Prp(x,y, z, ), is 

JL = x + y tanY - r - L secY - II-5 

P 

is the distance from P (x,y, z) to a point on the flat face of the 
detector, P^(x' = y%z') and J l is the distance travelled by the 

neutron from P (x' = 0, y' ,z' ) to the point of the nuclear reaction 

leading to detection. The light is considered to travel the normal dis- 
tance from the detection point to the phototube cathode. 

II. TIME COORDINATES 

In this study, the arbitrary convention is adopted that the zero 
of time for a given incident proton is at the instant it would have 
passed through the plane x = 0 if the target were absent. Apparatus delays, 

including those of the phototube, are ignored. These conventions do not 

affect the generality of the results. 



13 


If n is the conversion gain of the time analyzer in nanoseconds 
per channel, then the measured flight time t for an event is recorded 
in channel k if 


n(k - 1/2) < t m <; n(k + l/2) 


Il-k 


t is given hy 


t 


m 


t + t (h) + t. = t + t + t + t (h) + t. 
w j p n c w v/ j 


II-5 


where 


t is the total physical flight time, 

t (h) is the pulse-height dependent time slewing introduced by 
the detector, 

tj is the time jitter of the apparatus, 

t is the interval from the time zero defined above to the (p,n) 
interaction point in the target, 

t is the flight time between the birth and detection points for 
the neutron, 


t is the time required for light to travel from the neutron 


detection point to the cathode of the phototube. 


CHAPTER III 


MAXIMUM FLIGHT-TIME UNCERTAINTIES 

The time -of -flight of a neutron is related to the energy and the 
flight path in the nonrelativistic approximation by 


t (nsec . ) 


72.3 L(meters) 

1 7p > 

4 (MeV.) 


III-l 


and it follows directly that the relative resolution time for a given AL 
is simply (At/t) = (AL/l), where AL is introduced by the flight-path 
uncertainty between the birth and detection points within the target and 
detector. We shall calculate the maximum relative flight-time uncertain- 
ties At/t' in terms of the target and detector dimensions assuming that all 
dimensions have negligible effect except the one being evaluated; so that 
in the practical case the various effects would have to be combined. The 
convention will be adopted in which At is the difference between the 
maximum and minimum flight times while t = L/v^, where L is the distance 
from the center of the target to the center of the detector face and v^ 
is the neutron velocity. 


I. DETECTOR EFFECTS 

For the detector geometry shown in inset A of Figure 3, in which 
the radius produces the dominant flight-path variation, the maximum and 
minimum neutron flight times are given by 


Ik 
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t = — (l 2 + kL) 

max v \ T)/ 

n 


1/2 


III-2 


and 


t . 
mm 
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where is the detector radius. Expanding III-2 in a binomial series and 
subtracting III-3 yields the flight-time uncertainty At. Division by t 
gives the maximum flight-time resolution as 



III-4 


In typical applications, Rp/L < 0.1, and all terms beyond the first may be 
ignored since their total relative contribution to the resolution time is 
less than 10" 5 . 

Inset B of Figure 3 shows the geometry for a "thick" detector. 
Considering the effects of the detector thickness D, the maximum variation 
in the flight time is 



n 


KD 

c 


m-5 


17 


The second term on the right-hand side of the equation gives the correction 
for the light transit time through the scintillator medium under the 
simplifying assumptions given in the introduction. K is the index of 
refraction of the scintillator and c is the velocity of light in vacuum. 
Dividing both sides of III-5 by the nominal flight time t we obtain the 
maximum flight time resolution for the axial dimension of the detector as 
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Figure 4 shows the results of plotting the time resolution as a 
function of R^/L and D/L. Equation III -6 is plotted for two values of v 
corresponding to 2- and 100-MeV. neutrons. These data show the expected 
dominant effect of the detector thickness. 

It is interesting to note the relationship between Rp and D for a 
detector of constant volume and minimum relative flight -time uncertainty. 
The spectrometer efficiency r \ ^ is given by the product of the detector 
efficiency, e^, and the solid angle subtended at the target by the 
detector, is defined here as the ratio of the number of neutrons 
detected to the total number of neutrons incident on the face of the 
detector. Ignoring attenuation in the scintillator, £ (v^ )D, where 

S(v ) is the macroscopic efficiency cross section for detection in the 
sc int illator . Then, 



III-7 



1 


10- 5 


■e k. Flight -time 3 
s Rjj/L and D/L. 
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where is the area of the detector face. In this approximation 77 ^ is 
directly proportional to the detector volume. The expression for the 
maximum resolving time of detector in which both R and D are considered is 


At D 

■ _ 

t — L 


(a. 
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which is correct to the second order in R ^/L and D/L. Setting = VOL/mD 
and taking the derivative of III -8 with respect to D, we get the thickness 
for which (At/t) is a minimum, under the constraint of constant volume 
(given spectrometer efficiency), or 


D = 



III-9 


It follows from III-9 that for minimum At/t, R_ and D are related by 


Rj = 2DL 



III-10 


Equation III-10 is the relationship for which the radial and axial 
dimension contribute equally to At/t. Equations III-7 through III-9 yield 
the smallest maximum flight -time resolution as a function of the spectrom- 
eter efficiency to be 
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At 

t 


CVJ 


/ *LS(v n ) ' 

? « ) 
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It should be noted, however, that the relative dimensions specified by 
these equations will lead to poor angular resolution and/or inefficient 
light collection as one increases r^. In practical applications is 
therefore smaller than the value suggested by Equation III-10, so the 
thickness generally introduces the significant detector contribution to 
the geometric time resolution. 


II. TARGET EFFECTS 

The determination of the maximum flight-time resolution arising 
from the target dimensions requires a more detailed analysis. The flight- 
time uncertainties are a function of the rotation angles of the target and 
detector and the proton transit time through the target. 

Inset C of Figure 3 shows the thick target geometry. The protons 
are incident on the target along the x-axis in a line beam. The target 
thickness along the beam axis is S = T secY where Y is the rotation angle 
of the target and T is the target thickness. The extreme flight-times 
which arise for a given target -detector configuration depend on the 
particle velocity ratio v n /Vp and the detector rotation angle 9 . For 
forward scattering angles and v^ « v , the longest and shortest flight- 
times are obtained at the target edges as shown in the figure. However, 
at given detector angles and increased neutron velocity, the maximum and 
minimum flight-times may reverse because the proton transit time through 
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the target now comprises a significant fraction of the total flight time 
along a given path. In all cases when 6 > 90°, corresponding to backward 
scattering, the maximum flight time is the sum of the proton transit time 
through the target and the neutron flight time along the longest path from 
the target to the detector; this value is always greater than the neutron 
flight time along the shortest path to the detector measured from the 
target edge. 

Tha maximum and minimum flight times measured from the target edges 
are given with the aid of the cosine law as 

% = (L 2 + &fk + SL cose) 1 / 2 - g|- III-12 

n p 


and 


t x = — (L 2 + & /k - SL cose) 1 / 2 + , III-13 

n P 

where 6 is the detector rotation angle. S/v is the proton transit time 
through the target when v is the proton velocity. Expanding the square- 
root terms in III-12 and III-13 and following procedures outlined above^ 
we obtain the result that 

(At/t) g s 1 1 (cose - v n /v )| . III-14 


Equation III-14 is accurate through second order in S/L. 
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Depending on the detector rotation angle and the particle velocity 
ratio, given target -detector configurations arise where the maximum and 
minimum flight times given by Equations III-12 and III -13 are nearly equal 
For these cases, the minimum flight time occurs for interaction points 
within the target rather than at either target edge. Then, the total flight 
time is given by 


t(x) = — + (l 2 + x 2 - 2xL cos©) 1 / 2 HI-15 

P V n 

where t(x) is the flight time measured from any point on x. Taking the 
derivative of III-15 with respect to x and solving for x/L when dt(x)/dx = 0, 
one obtains 


v 

(f) — sin0 + cose ; HI-16 

' Ly t=min (v 2 -v 2 ) 1 / 2 
P n ' 

the value of (x/l) for which t(x) is a minimum provided v ^ v n * When 

v^ = v , there is a minimum at 9 = 0 independent of x. 

* \ 

Figure 5 shows the curves of (x/L), . versus 9 for 2- and 100-MeV. 

t=mm 

neutrons. These data give the range of values of e for which Equation 
III-14 is an invalid approximation for (At/t). For the range 77°-88°, 
there is a minimum within the target when S/L = 0.1. Over this range of 
angles where t (x ) < min (t x , t 2 ) Equation III-14 should not be used. The 
minimum flight time is then determined from Equation III-15 using the 
values of (x/L). and 6 obtained from Figure 5 or Equation III-16. The 


t - min 
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relative resolution time (At/t) is then calculated using the relation 
[max(t x ,t 2 ) - t(x)]/t. For the case at 0 given by (x/L) t _ m ^ n = 0, the 
flight-time resolution is simply 



Figure 6 shows the geometric induced flight-time resolution plotted 
as a function of Q for S/L =0.1 and 0.2 for 2- and 100-MeV. neutrons. 

The solid portion of each curve was obtained using Equation III-14. The 
dotted portion was obtained using Equation III-15; the range of angles 
being determined from Equation III-16 and Figure 5 for x/L just equal to 
the value at the target edges. 

The calculation of the effects on the flight-time resolution arising 
from the target radius follows similar arguments. Inset D of Figure 3 
shows the thin target configuration in which the radial dimension introduces 
the dominant flight -path uncertainty. In this case, the proton beam is 
taken as uniformly distributed over the target radius. 

Ignoring multiple scattering of the protons in the target, the 
"effective" target radius corresponds to the radial dimension of the beam, 
which in most applications is smaller than the target radius. If the 
cartesian components of the beam are studied separately, we observe that 
the z -component of the beam produces a AL which is independent of the 
detector rotation angle, 0. The calculation of the maximum timing un- 
certainties from the z -component yields a solution similar to 
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Equation III-4 obtained for the detector radius. The uncertainties are 
of second order in R_/L. The y-component of the beam, on the other hand, 
produces a AL which is dependent on 9 and Y and contains first order terms 
in B^/L. On the generally valid assumption that the beam radius is 
smaller than the detector radius, the maximum flight-time uncertainties 
can be calculated in the x-y plane with no practical loss in generality. 
For those cases where the beam radius is comparable to the detector 
radius, the contribution from the z -component can be added in quadrature 
with that obtained for the y-component. 

When the maximum and minimum flight times occur at the target 
edges, we obtain the expressions for the extremum flight times from the 
cosine law as 


and 


1 p/p Rrp sinY 

t 2 = 7 - [L 2 + B 2 - ZR^ sin (0 - Y)] 1 ^ - 

n p 


HI-18 


\ = “■ [L 2 + R^i + 2BipL> sin(e - Y )] l/2 + ^ 
n 


sinY 


II 1-19 


where B^ is the target radius. Expanding the square -root terms and 
dividing the difference | t g - | by the nominal flight time t yields 


\ 3 . 

2R t 

^'e t 

L 


n 


III-20 
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which is accurate through second order in R^/L. Equation III-20 is an 
approximation to the flight -time resolution which may he used provided 
the target edges yield the maximum and minimum flight times. For the 
other case, the actual minimum flight time exists for a point within the 
target . 

For neutrons horn at a point 77 in the target, the flight time is 
given hy 


t(r?) = ~ [L 2 + rf - 2rjL sin(e - Y )] 1 ^ 2 - % ^ inY . III-21 

V n V p 

Taking the derivative of III-21 with respect to 77 and solving for ( 77 / L ) 
when [dt(?])/d? 7 ] = 0 , we obtain the value of (tj/l) . for which t ( 77 ) is 
a minimum. The result is 



t=min 


v sinY cos (0 - Y ) 

775- + sin ( 0 - ¥ ) • 

(v 2 - sin 2 Y v 2 ) ' 
p n' 
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The equations which describe the target radius dependent time 
resolution are a function of the rotation angles 9 and Y as well as the 
particle velocity ratio. A typical case is discussed here. Figure 7 
shows the results of plotting ( 77 / L as a function of 6 for 2-MeV. 
neutrons. The target is rotated at 30° to the incident beam. Assuming 
B^/l = 0.1, we observe from Figure 7 that for the range of angles from 
28 ° to 39°, the minimum is within the target and Equation III-20 should 
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not be used to evaluate At/t. Instead the minimum flight time is evaluated 
using Equation III-21 when t(rj) < min(t 1 ,t 2 ) using the values of 77 /L and 0 
from Figure 7 for the given target rotation angle Y. If t ( 77 ) £ min^ ,,t 2 ) 
Equation III-21 is used to compute At/t where At/t = [max^^t ) - t(rj)]/t. 
For the case at 0 when ( 77 / L ) += Trn - ^ = 0 , the relative timing uncertainty is 
given by 



L 


sin (0 


v 

Y ) — sinY 

v 

P 
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The results of these calculations are illustrated in Figure 8 . The solid 
portion of the curve is obtained from Equation III-20 and the dashed 
portion using Equation III-21; the values for 0 and ij/l being determined 
from Figure 7* 






CHAPTER IV 


THE TIME CHANNEL RESPONSE FUNCTION 

The simplified spectrometer defined in Chapter I is analyzed in terms 

of the time-channel response function S^v jv n , e )> the probability per MeV. 

of neutron energy per incident proton that a neutron of speed produced 

in the target is detected in channel k of the multichannel-time analyzer. 

The time- channel response function is derived from the convolution of the 

flight-time distribution function and the time slewing response function 

and the integral over the analyzer sensitivity function. Each of these 

functions will be discussed in detail. 

S k includes the "response" of the target through its differential 

neutron production cross sections (v ;v ,Q ). The sum of S, over all time 

P P H 

channels is related to the counter efficiency by 


) S (v ;v * 0) as Z-Ay ) 
L> k v p n' 7 D n' 


A. 


D 


2 (v ;v , e)S 


P P n 


eff 
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where S is the target thickness along the beam direction taking into 
account the beam attenuation. 

I. THE PHYSICAL FLIGHT-TIME DISTRIBUTION OF DETECTED 
NEUTRONS AVERAGED OVER THE TARGET 

The flight -time distribution G is calculated here for neutrons 
which are born in the target and interact in the detector in terms of the 
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physical flight time of Equation II-5 . The distribution function is 
obtained by defining the geometric response function for neutron inter- 
action events in the detector as a function of the neutron flight time 
from a point in the target and then summing over the target volume . The 
flight-time density function of detected neutrons is the probability per 

incident proton of speed v per MeV. of neutron energy per nsec, that a 

]? 

secondary neutron of speed v^ has an apparent flight time t, and is 
given by 


G(v ;v n ,t) dt = dt JJ]“ N(y,z) B(A ,fi) ,x,v n ) dxdydz IV-2 


where dxdydz is an element of target volume and 


N(y,z) dydz = the probability that a beam particle has a 

displacement from the beam axis within dy at 
y and dz at 


B(j& dx 


the probability per unit solid angle per MeV. of 
neutron energy that the proton interacts in the 
target within dx at a penetration = x + y x 
tanY - 7*! sec Y to produce a neutron which is 
emitted in the direction Q, and 


Eg(t yl ,x,v ) dt = the detector geometric response function: The 

p n probability that a neutron of velocity v^ born 

at x,y in the target and emitted in a unit solid 
angle toward the detector is detected within dt 
at t. 


The Detector Geometric Response Function 

The detector geometric response function is given by 


E e (t,A p ,x,v n ) dt = dt J I(t,X ,x,v n ,CJ) dH 

a d 


IV-3 
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where the integration is performed over the area of the detector face, A^. 
We assume for all cases that the y' and z' dimension of the detector front 
surface are greater than those of the beam intersection on the target. 

<30 = [(!“’ • n)/X^.]<3A^ is the element of solid angle subtended by an 
element of detector front surface at the interaction point in the target 
P T (x,y,z). is the distance from P T to face of the detector and i' is 
the unit vector along the detector axis x' . The integrand, the probability 
that a neutron of velocity v n born in the target at x,y and emitted in 
the direction Q is detected within dt at t, is given by 


-fcpj x, v n ,Q ) dt 


( exp (" S TN f ^) Z(v n ) exp [' E D ff V^] d V t} 


provided the neutron 
interacts at a detector 
penetration distance 
An(t) along the direc- 
tion Q 


IV-4 


1 ° 


for j£j)(t) inconsistent 
with neutron interactions 
in the detector 


where 

eff 

E^ = effective macroscopic absorption cross section for 
neutrons in the target, 

E(v n ) = macroscopic "efficiency" cross section for neutron 
detection, 

cf f 

Ep = effective macroscopic absorption cross section for 
J neutrons in the detector. 


v n = neutron velocity, 

Xpp = pathlength traveled by the neutron in the target, 


3 ^ 


i (t) = pathlength traveled by the neutron in the detector 
expressed in terms of the apparent flight time, 

d£ (t) = differential distance traveled by the neutron in time 
D dt. 

The shape of l(t, x, v^, Q) in time is just an exponential truncated at 
the time limits corresponding to the time at which the neutron would pass 
through the front and rear surfaces of the detector. The combination of 
these exponentially shaped distributions for various points on the 
detector surface complicates the original simple shape. 

Once l has been written as an explicit function of t, Equation 
IV-3 gives the distribution of neutron flight times summed over the 
detector volume; the volume element in the detector having been defined 
in terms of dA^ and v^dt . 

The distance traveled by the neutron from its birth point in the 

target P^(x,y, z) to the point of detection P^ is t^v^. With the aid of 

Equation II-5 with t. = t (h) = 0, we can write 

J w 

tv = l + i (t ) = (t - t - t)v IV-5 

nn f D' /s c p'n 

where is the distance traveled by the neutron from P^ to a point on 
the face of the detector and j£^(t) is the distance from the detector 
face to the point of detection P^. It has been assumed that the light 
from the neutron interactions is normally incident on the photocathode. 
Then the light transit time from the point of neutron interaction to 
the photocathode is 


t 

c 


- 1 c d - V 1 ' 


■ CJ)] 


> 


I V-6 
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where (i/ . Q) is the angle at which the neutron strikes the face of the 
detector. Substituting Equation IV- 6 into Equation IV-5 and solving for 
yields the result 




[t - t p - (KD/c)]v n - £ f 

A 


lv -7 


and 


dJt D (t) = 


v dt 
n 

A 
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where A = [1 - (Kv^/c)(i^ • Q)]. As the denominator goes to zero, the 
only problems introduced are numerical and I is nonzero only over a small 
time hand. 

The proton transit time through the target, t , is obtained 

P 

assuming the time reference is at x = 0 with the target imagined to be 
absent. Then if dE/dx is considered constant through the target, 


t 

P 



i 2 

, £ 
4E v 
P P 
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where E is the proton energy (MeV.) and -(dE/dx) is the incremental 
P 

energy loss for protons in the target material. -(dE/dx) is a positive 

number. The derivation of t is given in Appendix A, while Appendix B 

p - 

contains the detailed analysis leading to expressions for (i' • Q), i , 
and A . 
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Integrals Over the Target Volume 

To obtain the complete flight -time distribution given by Equation 
IV-2, we must now integrate over the target volume taking into consider- 
ation all possibilities for proton interactions. In that equation. 


B(X p ,Q) dx = 


£ p (v p ;v n ,n) exp(-E^ f j& ) dx j x, y within target 


x, y outside target 


r p ( VV n) 


the differential cross-section per MeV. per unit 
solid angle in the direction Q for the production 
of neutrons of speed v^, 

the effective cross section for absorption of 
protons in the target. 


Since the chance of emission of the neutron is taken to be independent of 

the azimuthal angle, E (v ;v ,Q) reduces to E (v :v ,©), where © is the 

p p n' ' p' p n' ■” 

angle of neutron emission relative to the beam axis. If E (v ;v ,©) does 

p p n' 

not vary too rapidly with © or v , it may also be considered constant in 
the integral for G. [Experiments are normally designed to make E p (Vp jv^,© ) 
sufficiently constant over the detector face so that E^C©) — £ p (0)]« 

It has been shown in Chapter III that the z -component of the beam 


typically introduces negligible flight-path variations. Based on these 
results, Eg(t,i,p,x,v n ) and B(i.p,Q) were considered independent of z, and 
since N(y, z) is the only term in which z is contained explicitly, it is 


integrated immediately over a] J values of z with the result that 
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dy J* M(y,z) dz = N(y) dy . IV-10 

N(y) is approximated by a gaussian distribution with mean value equal to 
zero, given by 


S (y) = -i_ e (-i /2)krfr? 


IV-11 


where 

ct = the standard deviation of the proton beam distribution about 
the x-z plane. 

II. THE EFFECTS OF INSTRUMENTAL TIME JITTER 

There is a spread in the duration of the time intervals between 
the production of photoelectrons resulting from the neutron interaction 
in the detector and the time when the amplitude of the output pulse 
triggers the electronic circuitry. This time spread depends on several 
factors including the energy left behind by the particle in the scintil- 
lator, the conversion efficiency of the scintillator, and the statistical 
variations in the cathode photoelectron emission and in the electron 
multiplier tube. The magnitude of the timing jitter depends on the 
fluctuations in all of these causes as well as the amplitude of the pulses 
introduced by the light output in the detection event. The amplitude 
variations are treated in Section IV of this chapter and the time jitter 
is handled separately under the assumption that it is independent of 
amplitude variations in the output pulse. 
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Included in the jitter distribution are the fluctuations arising 
from the distribution of the protons in time and the time variations 
introduced by the zero-time detector. These uncertainties, along with 
those described above, are here lumped into one constant source of 
instrumental time jitter. 

The jitter manifests itself in a broadening of the apparent flight- 
time distribution. The causes of jitter arise randomly in time so that 
the time-jitter distribution can be approximated by a normal frequency 
function with zero mean jitter time and constant standard deviation o' ., 

O 

given by 

J(t.) = -i— . Iy-12 

The jitter-smeared apparent flight time t^ for an event is 

t . = t + t . . IV-13 

i J 

Then, the jitter- smeared flight -time distribution function is obtained by 
the convolution of Equation IV-12 with IV-2: 


: i = dt i S J(t i " 


17-14 


Completing the integration of Equation IV-14 proceeds as follows. 
Equation IV-7 is inserted into Equation IV-4 and all terms containing the 
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physical flight time t are factored out. These terms, along with Equation 
XV-12, are placed inside the integral sign in Equation IV-14. Completing 
the square of the terms in t and integrating over all values of t yields 
the following result 


at i C( W V - dt i T n S( V S p ( VV e) 0 B <V S) 


eff . 
e ™ T /A 


x exp[(^ ff /Aj u f + KDv/c + t p T n )] 
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X exp 



W2 


)] 

l"erf j 

1 

L ' 



erf 


it . 
mxn 


a 


j 



dDdxdy , 


where 


A = 1 - (Kv n /c)(i' • Q) 


IV-15A 


and 



2 „ef f / . 


IV-15B 


t . and t are the minimum and maximum physical flight times t along 
mm max 

the direction Q to the front surface and point of exit of the detector. 


t . and t are given by 
mm max 


t . 

mm 


tp + 4 f /v n + KD/c 


and 
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W = *p + V v n + (D/v n V 0 -' ■ S > 

from Equations 17-5 through IV-7* The error function used here is defined 

by 



,-x?/ 2 a a 


dx = erf — 
a 




erf — 

a 


17-18 


In practical calculations, the quantity C(v ;v n ,t^)/E (v ;v n , 0) is 
evaluated since the purpose is to calculate Ep(v ;v n , 0) from experimental 
results. 


III. RELATIONS LEADING TO THE DETECTOR EFFICIENCY 

The integral of Equation IV-14 over the jitter-smeared flight time 
t^ should yield the efficiency of the detector times the number 5 of 
neutrons which strike the face of the detector, or 


It. dt i It at G < v p ; V t > J ( t i - - e d s 
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» = yVp, v e)T eff r aH . 

JD 


where 


hi 


T eff is the effective target thickness given by 


T 


eff 
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if the neutron attenuation in the target is ignored. In this section, 
we shall perform the integration indicated in Equation IV-19 in order to 
prove the validity of the relationship. 

Changing the order of integration in Equation IV-19 and performing 
the inner integration yields the result 


Jt Q (W t} dt " e D 5 * IV-21 

It is possible to perform the integration in this manner since j(t ) is 

j 

normalized to unity. The convolution of any function with one of unit 
area yields a function with the same area as that under the original 
function. 

In the assumption that no neutron attenuation occurs in the target. 
Equation IV-4 reduces to 

dt I' (t,A p ,x,v n ,5) = 2(v n ) expC-Ejy^t)] <U D (t) . 17-22 

Since this is the only equation in which t is contained explicitly, it 

may be integrated over t immediately. Substituting ) and dAp(t) 

from Equations IV-7 and IV-7A into IV-22 and integrating from t to t 

mm max 
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We obtain the result that 
t 



where t . and t are given by Equations IV-l6 and X7-17> respectively, 
mm max 

and D/(i' ♦ fi) is the distance traveled along the direction 0 by the 
neutron in passage through the detector. 

Now, with the aid of Equation 17-2, we can write 

€ T? - Xy Xx(y) Xa,, e D ( V X ’V Sj N(y) BU p ,5) dfldxdy . IV-24 

e^(£ ,x,v ,fl) is a function of the coordinates of the neutron birth 

D p n 7 ' 

point in the target and the direction of emission of the neutron. If 
the direction of emission of the neutron is such that it passes through 
the detector front surface, then Equation 17-4 is nonzero. Further, the 
spectrometer efficiency, given by the product of the detector efficiency 
and the element of solid angle subtended by the detector, does not 
depend strongly on the neutron birth coordinates. Then Equation 17-24 
may be integrated as follows 
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e D (v n ,Q)dQ 


target 


J y J x ( y ) e D (t p ,x,v n ,n) W(y) BU p ,Q) dQdxdy 

IV -25 

Jy T x(y) N(y) BU p ,Q) dx dy 


where ^ J e^(v n ,Q) target is the average spectrometer efficiency 
weighted by the interaction density in the target. In this analysis, 
is considered essentially constant over the detector area, so that 


e D (r n ) " 


Jaj, e D ( V S > «) 1 


target 


iv-26 


AQ 


e i)(Vn ) is ^ ust tlie avera S e value of the detector efficiency and = 


Ja^ dn is the element of solid angle measured relative to the origin of 
the target coordinate system. Now, we can write 


e D (v n ) An J f N(y) bU ,e) dxdy l iv-27 

y 'x(y) * 

since ^p( V p-’ v n >^) ^p^p^n' 0 ^’ T ' ie integration in Equation IV-24 has 

now been greatly simplified. 


bb 


The limits of integration for the x and y integrals are given by 


T 1 secY - y tanY £ x ^ r 2 secY - y tanY 


IV-28 


and 


-2.5a ^ y ^ 2.5a 


IV-29 


Introducing IV-9 and IV-11 into IV-27 and completing the integration 
over the specified limits, we obtain 


e (v ) 
' n' 


£ (v 
P v P 


;v ) e ) 

n ’ ' 
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Then, the integral of Equation IV-19 over the physical flight time yields 
the average value of the detector efficiency weighted over the target. 

IV. THE TIME-SLEWING RESPONSE FUNCTIONS 

The large amplitude range in the light output from the neutron 
interaction causes a portion of the pulse amplitude spread to be trans- 
lated into a timing spread which is called slewing or walk. The time- 
slewing response function W[t^(h), v^] dt^(h) is defined as the probability 
that a neutron of velocity v^ produces an amplitude- dependent slewing 
time within dt (h) at t w (h). The response function is normalized such 
that the integral over all values of t^(h ) is unity, 

f W[t (h),v ] dt (h) = 1 . 

J t (h) w n w 
w' * 


IV-31 
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In order to determine the slewing response function, the pulse- 
height spectrum resulting from the neutron interactions in the detector 
and the time-slewing curve must be known. The pulse-height spectrum is 
calculated using a first collision approximation to obtain the differ- 
ential macroscopic cross section. Conversion of these data to the 
corresponding pulse -height spectrum is accomplished using the fast -light 
output curves for charged -particle interactions in the detector. The 
slewing curve which gives the time shift in terms of the ratio of the 
pulse height to the amplitude of the output pulse which just triggers 
the time analyzer, was taken directly from the data of Peelle etal. 24 

Derivation of the Pulse-Height Spectrum 

Neutron detection is achieved in organic scintillators by the 
elastic scattering on hydrogen nuclei and the nonelastic interactions 
with carbon nuclei. Both give energetic charged particles which produce 
light in the scintillator. 

The cross sections for these reactions will be derived using the 
method and notation of Sc hut tier, 8 5 and the resulting energy spectra 
converted to pulse-height spectra in light-output units using the fast 
light-output curves of both Schuttler and Birks. 26 

24 R. W. Peelle et al.. Differential Cross Sections for the 
Production of Protons in the Reactions of 160-MeV. Protons on Complex 
Nuclei, ORNL-5887, September 1966. 

25 R. S. Schuttler, Efficiency of Organic Scintillators to Fast 
Neutrons, ORNL-3888, July 1966. ~ 

26 J. B. Birks, Proc. Phys. Soc. A 6k, 74 ( 1952 ). 
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In this analysis, we are assuming that the charged particles from 
the hydrogen reactions are totally absorbed in the scintillator. For 
example, for 2.5^--cm. thick detectors, the assumption is valid for 
neutrons up to lk.5 MeV. where the maximum range of the charged particle 
is 0.2J cm. For 50-MeV. neutrons, however, the range of the charged 
particle is 2.20 cm. which is comparable with the detector thickness, 
and the resulting spectrum is in error. 

Hydrogen scattering. In the energy region below 10 MeV., (n, p) 
scattering is isotropic in the center-of-mass system and all proton 
energies from down to zero are equally possible after a single scat- 
tering. At neutron energies greater than 10 MeV., the scattering becomes 
increasingly anisotropic . 

The energy differential macroscopic cross section for elastic 
scattering of a neutron of energy in hydrogen is given by 

2 h (e n ;E p ) = n h (e n cos0)_1 iv “ 52 

where <j(Ejj, 0) is the laboratory microscopic differential cross section for 
scattering at angle 6 and N„ is the atomic density of hydrogen in the 
scintillator. a(&) is approximated using semiempirical formulae due to 
Gammel. 27 The relationship between E^ and 6 is obtained from classical 
kinematics as 


27 J. L. Gammel, Fast Neutron Physics, Part II, (Interscience 
Publishers, New York, 1963), p. 2209* 


^7 


E p = I ^ ^ ‘ cos20 ) = ^ sin 2 e • 17-33 

The energy differential cross section is obtained for values of 6 

in the range 0 ° < 0 ^ 90° • Equation TV-33 yields the value of E which 

P 

when replaced in IV-32 gives the value of ^(E^E^). 

Figure 9, taken directly from the work of Schuttler, 28 shows the 
proton-energy differential macroscopic cross section for elastic scattering 
on hydrogen plotted as a function of the proton energy for several values 
of the neutron energy. 

Reactions in carbon. At neutron energies greater than 10 MeV., 
the probability for nonelastic scattering with carbon becomes important. 

The most important ground state reactions which compose the nonelastic 
cross sections and which give charged particles are 

1. C(n,n*)C* leading to C* -* 3 ot; Q 5: -7-8 MeV., 

2. 12 C(n,p) l3 B; Q = -12.6 MeV. and l3 C(n,np B; Q = -l8A MeV., 

3. 12 C(n,ctf Be; Q = -5.7 MeV. 

The cross sections for the reactions described in 1 and 2 above have been 
calculated by Schuttler 29 using an adaptation of the Monte Carlo code of 
Bertini 30 which is based on the intranuclear cascade model. Schuttler's 
analysis did not produce ground state reactions (n,a) leading to 9 Be. 


28 Schuttler, op. cit., p. 21. 
29 Schuttler, op. cit., p. 28. 


30 H. Bertini, Phys . Rev . 131, 
Calculations on Intranuclear Cascades 


1801 (1963 ); also Monte Carlo 
, OENL-3383 (1963). 
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We shall adopt his method for deriving the differential cross section 
for this reaction. 

It is assumed that the reaction products from the (n,a) reaction 
are isotropic in the center-of-mass system. The differential macroscopic 
cross section for the production of alpha particles is then given by 


N n cr (E - E . T" 1 ; 
C a amax amin' * 


J C(n,a) \ 


0 


E . < E < E 

a nun a a max 

iv-34 

otherwise 


where <7 is the total cross section for this reaction, N_ is the atomic 
u, C 

density of carbon in the scintillator, and E arrax and E anin are the maximum 

and minimum possible a -particle energy values. E and E . were 

amax amin 

calculated using the momentum relations assuming that only the neutron is 
relativistic. The cross section is then taken as being uniform in the 
corresponding energy region. This approximation is reasonable since 
Oq, ^ 0.1 a c (nonelastic) over the allowed energy region for which the 
reaction is possible and the contribution to the total pulse-height 
spectrum is small compared to the hydrogen and other carbon reactions. 

Figure 10 shows the relative differential cross section for the 
cascade reaction plotted by Schuttler 31 as a function of the electronic 
bias (in light units) for a 2.5^-cm. thick scintillator. The detailed 
analysis leading to this cross section is given in ref. 20. For the 
purpose of this analysis, we derive a standard energy spectrum based on 


31 Schuttler, op. cit ., p. 36. 


2-01-058-1018R 


a 0.05 



Eg (MeV e.q. 


Figure 10 . The effective relative differential cross sections 
for a 2.54- cm^thick NE-213 scintillator as a function of the bias 
energy (electron equivalent). These data are from R. J. Schuttler, 
Efficiency of Organic Scintillators for Fast Neutrons, ORNL-3886> 
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the data of Figure 10 which, when combined with the appropriate effective 
eff eff 

cross section for the cascade reaction, yields a rough 

estimate of the pulse-height spectrum for the cascade reactions. 

To obtain the pulse-height spectrum for the interactions of neutrons 

of any energy, all that is required is to multiply the ordinate of the 

"standard" spectrum shown in Figure 11 by the appropriate cross section 

eff 

from Table I, where Table I gives the values for ct for the reaction 
as a function of the neutron energy for two detector thicknesses. 

The total spectrum . The final step in the derivation of the total 
pulse-height spectrum is to convert the proton and alpha-particle spectra 
given by Equations IV-32 and TV-34, respectively, to light-output units. 

The fast light -output curves for protons and alpha particles 
plotted as a function of neutron energy are given in Figure 12. These 
data are taken from Schuttler 32 and Birks. 33 The light output is measured 
in terms of an electron which produces the same amount of light as the 
charged particle being measured. Converting the energy scale for the 
charged-particle reactions to electron-equivalent light output is accom- 
plished directly from these curves. The cross section at each value of F, 
however, must be divided by the slope of the light curve where the slope 
of the curve is given by Birks 34 for the fast light output as 

32 Schuttler, op. cit., p. 53* 

33 Birks, op. cit., p. 75* 

34 Tbid. 
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E c (MeV ELECTRON EQUIVALENT) 


Figure 11. Standard energy spectrum for the carbon cascade 
spectrum. 
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TABLE I 

-P-P ±L 

EFFECTIVE CROSS SECTION AT E^ = l80 keV. e. q. 


\ (MeV.) 

eff , , x 
0 q (nib . ) 

For I) = 6.1 cm. 

For D = 2.5^ cm. 

10 + 

94 

65 

12 

128 

91 

1^4- 

187 

150 

15 

207 

168 

20 

319 

313 

25 

326 

320 

30 

334 

331 

40* 

310 

310 

* 

These data are 

from R. J. Schuttler, 

Efficiency of Organic 

Scintillators for Phst 

Neutrons, ORNL-3888, 

July 1966. Ep, is the bias 


setting of the detector discriminator. 180 keV. electron equivalent 
(e. q. ) corresponds to 1-MeV. neutron energy. 

^ Belov ~7.8 MeV., the cross section goes to zero. 

* 

For neutron energies greater than 40 MeV. the cross section is 
nominally 310 mb. for both detectors. 


LIGHT OUTPU! (MeV) 
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Figure 12. Fast light output for protons and alpha particles. 
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dF/dE = [1 + B(dE/dx)] -1 . IV -35 

B is a constant relating the energy loss of the charged particle to the 
light output (the value of B is taken as 0.0140 g. cm. -2 MeV." 1 in 
obtaining the curves in Figure 12) and dE/dx is the energy lost per unit 
pathlength by the charged particle in passage through the scintillator. 

When all the spectral data have been converted to electron- 
equivalent light -output units (cm. -1 MeV." 1 electron equivalent) all that 
remains is to add the proton recoil and alpha-particle spectra with the 
cascade spectrum. Figure 13 shows the resulting pulse-height spectra for 
the interactions of 2- and 1^.5 -MeV. neutrons in the NE-213 scintillator. 
Also shown is the bias setting for 1-MeV. neutrons. Plotted for the 
purpose of comparison is the 1^.5 -MeV. spectra determined experimentally 
by Love et al . 35 

Conversion to Time Coordinates 

The time-slewing response function is obtained by converting the 
total pulse-height spectrum to time coordinates using the time-slewing 
curve. In this analysis, we are using the time-slewing curve of Peelle 
and his coworkers 36 which is modified to be consistent with the defini- 
tions of the neutron flight time discussed in Chapter II. The time- 
slewing curve, shown in Figure 1^, plots the ratio of the pulse height 

35 T. A. Love et al., Absolute Efficiency Measurements of NE-213 
Organic Phosphors for Detecting 14. k- and 2.6-MeV. Neutrons, ORNL-3593> 
September 1966. — 

Ss R. W. Peelle et al.. Differential Cross Sections for the Produc- 
tion of Protons in the Reactions of 160-MeV. Protons on Complex Nuclei, ~ 
ORKL-3«&7, September I966. 



SPECTRUM 
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6 5 4 3 2 ( OH 

SLEWING TIME (nsec) 

Figure l^. The time -slewing curve for the neutron detector. 
These data are from R. W. Peelle et al ., ORNL Report No. 5887. For 
this analysis the data have been modified to be consistent with the 
definitions for the time and may appear reversed when compared with 
the original data. 
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to the time analyzer discriminator trigger level as a function of the 
time shift or slewing. 

The zero slewing time shown in the figure is taken arbitrarily. 

In an experiment, this zero value must be determined specifically for the 
detector and instrumentation being used. The "zero" flight time is 
initially established from measurements of the apparent flight time of 
the highest energy particle available. For experiments of the type being 
studied here, the proton beam is usually used. The time slewing as a 
function of the pulse amplitude for constant flight time is then measured 
and the zero time slewing is assigned for pulses of the same amplitude 
as those occurring during the determination of "zero" flight time. 

Converting the pulse -height spectrum to time coordinates is a 
straight-forward process. The pulse height (abscissa of Figure 13 ) is 
initially divided by the discriminator bias setting, taken here to be 
90 keV., and the corresponding time slewing is read directly from 
Figure 14. The ordinate of the pulse-height spectrum is then multiplied 
by the slope of the walk curve where the slope is defined in units of 
(MeV./nsec. ). 

Figure 15 shows the slewing response function for 2- and 14.5-MeV. 
neutrons. These data have been normalized to unity consistent with the 
definition given by Equation XV-31* The abscissa is in units of time and 
the ordinate in reciprocal time units. It should also be noted that both 
curves show a maximum slewing time of 4.3 nsec. This corresponds to the 
walk at a detector bias setting of 1 MeV. for neutrons (l80-keV. electron 
equivalent); the value of the bias used in all calculations. 
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V. THE TIME -ANALYZER RESPONSE FUNCTION 

The time -analyzer response functions are the probability that a 
neutron of measured flight time t (see Equation II-5 ) produces a count 
in channel k of the multichannel time analyzer. The analyzer response 
function is approximated by a rectangular distribution over the width of 
a time channel by 


W = { 


0 ; 


V 


n ( k - D < v s n ( k + 1) 

otherwise 


IV-36 


where n is the analyzer conversion gain in nsec ./channel. 

VI. INTEGRALS LEADING TO THE TIME -CHANNEL RESPONSE FUNCTION 

The probability per MeV. of neutron energy per incident proton that 

a neutron of speed v is detected in channel k of the multichannel time 

n 

analyzer is the time -channel response function and is given by the 
expression 


S, (v ;v ,Q) = f A. (t ) T(v ;v ,t ) dt IV-37 

k v p' n S, * * * 9 ' nr v p* n* nr m 

m 

where A^(t m ) are the time-analyzer response functions and T(v jv^, t^) dt^ 
is the probability per incident proton per MeV. of neutron energy that a 
neutron of speed v has a measured flight time within dt at t . 

T(v^;v n , is obtained from the convolution of the flight -time 
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distribution of detected neutrons averaged over the target volu me , 

C( V p-> v rr , *'i)-> an< ^ time-slewing response function, W[t (h),v ], given by 

"w n 


T <VW “ J ^V'Vh) »(»■ - W dt l 

i 


17-38 


where t w (h) = t ffl - t from Equations II-5 and 17-13. 

The summation of 3 k ( v p ^' v n ; 0 ) over all the time channels yields 
the average detector efficiency e^v^) times the number of neutrons which 
strike the detector given by Equation 17-30. With the aid of Equation 
17-3 6, we can write Equation 17-37 as 


h ( vv)-EJ. 


n[k+(l/2)] 

1 • T(v ;v ,t ) dt 
n[k-(l/2)] P n m m 
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Using the same arguments which led to Equation 17-21, it can be shown that 


r T(v ;v ,t ) dt = r C(v ;v ,t.) 

J. ' p n’ nr m J. ' p’ n’ i' 


dt . 
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m 


t . 

l 


Introducing the results from Section III for the integral of C(v :v , t. ) 

p n* ir 

over the jitter -smeared flight time, we obtain the result 
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CHAPTER V 


COMPUTATION OF THE TIME -CHANNEL RESPONSE FUNCTION 

The equations leading to the time channel response function dis- 
cussed in Chapter IV are sufficiently complex that a computer was required 
to facilitate the computations. Four computer programs were written, each 
designed to complete a separate phase of the analysis. These programs are 
discussed briefly in Figure 1 6 which shows, in block diagram, how they are 
combined to produce the desired response functions. In this chapter, the 
methods of analysis used in each program will be discussed along with some 
of the particular calculations which are necessary to understand the solu- 
tion of the problem. The specific details and descriptions of the programs 
and instructions for their use are given elsewhere. 37 

I. THE NUMERICAL CALCULATION OF THE JITTER -SMEARED PHYSICAL 
FLIGHT-TIME DISTRIBUTION OF DETECTED NEUTRONS 

The numerical calculation of the jitter-smeared physical flight - 
time distribution was accomplished by completing the integration shown in 
Equation IV-15 and evaluating the results at values of t^ for which the 
distribution is nonzero. The integration was performed using Gaussian 
quadratures with all computations being performed on a Control Data 
Corporation 1604 Computer. 


37 R. T. Santoro, A Fortran Program for Calculating the Time 
Channel Response Function in Neutron Time-of -Flight Spectroscopy; 
to be published. 
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Figure 16. Block diagram and description of the computer 
programs for calculating the time-channel response function. 
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Detailed discussions of the Gaussian quadrature procedure are given 
by Lanczos 38 and Kou. 39 This method of integration was chosen since it 
yields an accurate approximation to the integral and requires fewer com- 
putational steps than other methods. In this section, we shall discuss 
the more important aspects of the calculations and some of the criteria 
used to evaluate the results . 

Tabulation of Cross Sections 

The calculations made in this analysis were based on the use of a 
carbon target and an NE-213 scintillator. The assumption was made that 
the effective macroscopic cross section for neutron absorption in the 
target is zero. For thin targets, in which the neutron mean free path is 

much greater than the distance travelled by the neutron from the (p, n ) 

interaction point to the point of exit from the target, the assumption is 

quite reasonable. For example, for 2-MeV. neutrons the total cross section 

for carbon is ~1«7 b. and the corresponding mean free path is ~7*5 cm * 

The effective macroscopic cross section for neutrons in the 
detector was obtained from 

- Vh + Vc,ne ’ V ’ 1 

where N and N^ are the atomic densities of hydrogen and carbon in the 
detector, ^ is the microscopic total cross section for hydrogen, and 

II 

38 Cornelius Lanczos, Applied Analysis, (Prentice Hall, Inc., 
Englewood Cliffs, New Jersey, 1956), pp. 396-^10 . 

39 Shan S. Kuo, Numerical Methods and Computers, (Addison -Wesley, 
Reading, Massachusetts, 1965 ), pp. 239-252. 
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a c ne is the nonelastic microscopic cross section for carbon. Table II 
lists the values for and ng over the range of neutron energies from 
2 to 100 MeV. 

The macroscopic efficiency cross section £(v ) was taken from 

n 

Schuttler. 40 These data are based on first collision calculations of the 
detector efficiency for bias values of 180-keV. electron equivalent. This 
corresponds to 1 MeV. for neutrons. Schuttler' s data extend up to 50 MeV. 
For higher neutron energies, the efficiency cross section was inferred 
from the total hydrogen and carbon cross sections. Figure 17 is a curve 
of S(v n ) versus neutron energy. Also shown in the figure are macroscopic 
cross sections for hydrogen and carbon which contribute to the efficiency 
cross section. 

The macroscopic absorption cross section for protons in the target 

eff 

£,pp was estimated using the nonelastic cross section for neutrons in 
carbon from the data of Bertini. 41 was taken as 233 mb. This value 

compares favorably with the value given by Voos and Wilson of 225 mb. 42 

Reduction of the Integral to Quadratures 

The evaluation of Equation IV-15 requires the solution of a four- 
fold integral evaluated at the values of the jitter-smeared apparent flight- 
times t^ for which the distribution is nonzero. The range of flight times 
for which the distribution is nonzero is implied by the conditions of 

40 R. S. Schuttler, Efficiency of Organic Scintillators to Fast 
Neutrons, ORNL-3888, July 1966. 

41 H. Bertini, Phys. Rev. 131, 1801 (1963). 

42 R. G. P. Voos and R. Wilson, Proc. Roy. Soc. (London), A23 6, 

to (1956). s b — ^ 
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TABLE II 


VALUES OF cr„ 
11 

AND cr_ 

C, ne 

FOR 2 <; Epj <: 100 MeV. 


* 

a H 

CT C,ne 

MeV. 

barn 

barn 

2 

2.90 

1.72 

5 

1.63 

1.15 

10 

0.9k 

0.64- 

15 

0.65 

0.87 

20 

0.4-8 

1.04- 

25 

0.38 

0.92 

30 

0.31 

0.87 

35 

0.26 

0.82 

4-0 

0.22 

0.76 

50 

0.17 

0.64 

60 

0.13 

0.53 

70 

0.10 

0.43 

80 

0.09 

0.35 

90 

0.08 

0.30 

100 

0.07 

0.26 


Neutron Cross Sections, BNL-325> Second Edition, Brookhaven 
National Laboratory, February 1966. 

^ M. Kalos and H. Goldstein, Neutron Cross-Section Data for 
Carbon, NLA 12-16 (1956). 


EFFECTIVE MACROSCOPIC CROSS-SECTION (cm 
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Equation IV-4. The values of used in the evaluation of the integral 

were determined from the extremum flight paths for the neutron for the 

particular target -detector geometry being studied. The minimum and 

maximum flight times were defined by (t.) = (A ) /v - 2.5 a. and 

1 min 1 min n 3 

(t.) = (i„) /v + 2.5a. where (i,„) and (£„) are the extremum 

1 max 1 max n J 1 min 1 max 

flight paths for the neutron between the target and the detector and a • 

is standard deviation of the jitter distribution. The addition and 

subtraction of 2.5 a. insures that the values of t. are defined well into 

J i 

the region of the tails of the jitter smeared distribution. The time 

mesh was then taken in equal intervals between (t.) and (t.) 

1 min 1 max 

The integral over the face of the detector was evaluated by 
expressing the element of detector area in terms of the cylindrical 
coordinates on the face of the detector. Then 

dA D = pdpda V-2 

where p and dp are the radial components and da is the element of azimuthal 
angle. The integration over the face of the detector was broken down into 
two regions in order to take into account the effects of the detector 
edges. The contributions to the flight -time distribution were then added 
according to their respective weights and the integral was evaluated over 
a as follows 

«/2 ^ jt/2 *4) 

2 f f pdpda = 2 P da f pdp + f pdp V-3 

-”/ 2 0 Lo %( a ) J 
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where is the value of the radius for which the projected flight-path 
of the neutron along Q just intercepts the point where the rear surface 
and the cylindrical surface intersect. The factor of two in the 
integration is required since the integrals are evaluated over half of 
the detector. The integration over the target volume is straight- 
forward; the limits of integration for x and y being given previously 
by Equations IV-28 and IV-29. 

The most tedious problem faced in computing the results was the 
determination of the proper number of quadrature points required in the 
integration. This was accomplished in several ways. First; the results 
were normalized to the calculated area under the Gaussian beam distribu- 
tion; the area being evaluated using the same number of quadrature points 
as in the integral over y . This procedure reduced the sensitivity of the 
distribution to the number of quadrature points in y. Then; the area of 
the detector face was inspected for a given number of quadrature points 
in p and a. If the computed area was the same as the known area of the 
detector face; the computation was allowed to proceed. 

These steps; however; did not always yield the optimum approxima- 
tion to the solution of the equation. The use of very narrow jitter 

distributions; where <j . was small compared to At.; sometimes required 

J i 

a larger number of quadrature points and/or a tighter time mesh. Often, 
the best results were reached by trial and error and the correct number 
of quadrature points was determined when the value for the integral 
converged and increasing the number of points no longer affected the 


results . 
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Testing the Validity of the Results 

After inspecting the results to insure that the solution to the 
integral had converged, the resulting jitter-smeared flight -time 
distribution was integrated over all values of t^ to obtain the detector 
efficiency. The solution was assumed correct if the efficiency calcu- 
lated in this manner compared to the efficiency calculated from 
Equation IV-23. This equation was computed using the values for the 
maximum and minimum flight times used to obtain the distribution. 

II. NUMERICAL CALCULATION OF T(v ;v ,t ) 

p n' m 

The jitter-smeared physical flight-time distribution and the 
time-slewing response function were convoluted as shown in Equation 
IV -38 using trapezoidal quadratures in Program C0NV0L. In this calcu- 
lation, the time -slewing distribution was divided into equal time 
intervals such that the initial and final time -slewing values formed 
the upper and lower time limits. The physical flight-time distribution 
was then divided into intervals of the same size. Forming the intervals 
in this manner insured that the edges of the time-slewing distribution 
(see Figure 15 ) were properly represented in the quadrature. 


III. INTERPOLATION AND ENTABULATION 


In this analysis, we have made wide use of the TERP code of 
Putzulu et al. 4 3 This code enables the user to make rapid accurate 
numerical approximations to functions which are not normally provided 
in the computer libraiy. 

All of the cross-section curves and tables of the error function 
were approximated using this system. The system yields specific values 
for these approximated functions by parabolic interpolation within the 
tabulated values supplied to the system. 


43 D. T. Putzulu, P. Aebersold, and W. R. Burrus, Description of 
the TERP System, A Fortran II/ 63, IV/63, 360 System for Efficient* 
Interpolation, ORNL-TM-I7O6 (December b, 1966). 




CHAPTER VI 


RESULTS 

In this chapter, we shall present some typical results from the 
numerical solutions of the equations presented in the preceeding chapters. 
Some characteristic jitter smeared physical flight-time distributions of 
detected neutrons averaged over the target will be illustrated for the 
interactions in the detector of 2-, l4.5-j> and 100-MeV. neutrons and the 
effects on these distributions arising from the instrumental time jitter 
will be demonstrated. 

Comparisons will be made between the maximum flight -time uncer- 
tainties induced by the target and detector dimensions using the equations 
derived in Chapter III and the uncertainties predicted by the computed 
flight-time distributions for the same conditions. Finally, the effects 
on the flight-time distributions arising from the convolution with the 
time-slewing response function will be discussed and some typical time- 
channel response functions will be presented. 

I. COMPUTED PHYSICAL FLIGHT-TIME DISTRIBUTIONS 

The physical flight -time distributions for detected neutrons 
averaged over the target were obtained from the numerical integration of 
Equation IV-15. For the purpose of this analysis, the target -detector 
configuration and parameters summarized in Table III were chosen to 
demonstrate the results ; these parameters typifying a practical choice 
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TABLE III 

PARAMETERS USED IN COMPUTING THE PHYSICAL 


FLIGHT-TIME 

DISTRIBUTION 



Proton Beam 

Target Orientation 

E n = 160 MeV. 

Composition 

Carbon 

ct = 0.25 cm. 

Thickness (T) 

0.5 cm. 


Rotation Angle (y ) 

0 

0 


T + T 

0 


X 2 

Detector Parameters 


Normal Distance from x = y = 

z = 0 to the 

70 cm. 

Face of the Detector (l) 

Neutron Scattering Angle (9) 


o\ 

0 

0 

Thickness (d) 


6.01 cm. 

Radius (R^) 


6.0 cm. 

Index of Refraction (K) 


1.^3 



Based on measured values used in the parent experiment. 
See R. T. Santoro, The Space, Time, and Energy Distributions 
of the Proton Beam of the Harvard University Synchrocyclotron, 
ORNL-3722 (January 1965;. 



in the light of possible experimental arrangements for the target and 


detector . 

The distributions were approximated using four -point Gaussian 
quadratures in x, y, p , and a with a 40 -point mesh in t^ for secondary- 
neutron energies of 2-, 14.5-, and 100-MeV. The results of the computa- 
tion are shown in Figure l8 for the cases where a ■ = 0.43 and 0.86 nsec. 

J 

The abscissa gives the jitter smeared apparent physical flight time and 

the ordinate is the probability per incident proton per MeV . of neutron 

energy that a secondary neutron of speed has an apparent jitter 

smeared flight time per nsec, at t_^. The ordinate values have been 

divided by the effective target thickness, the detector solid angle, and 

E (v :v ,©). The resulting distributions have been normalized to the 
p' p n 7 ' 

computed area under the beam distribution. 

All of the distributions with the exception of the one computed 

for 2-MeV. neutrons with o. = 0.43 nsec, are essentially Gaussian in 

J 

nature with the chief contribution to the time spread arising from the 

instrumental jitter. For example, atl4.5 MeV, the mean flight time for 

the neutrons over the 70-cm. flight path is 14.15 nsec. The resolution 

time assuming only presence of the jitter is 7*1$> an( i l4.2fo when a. = 0.43 

J 

and 0.86 nsec., respectively. 

For the single case where E = 2.0 MeV. with a . = 0.4-3 nsec, the 

n 0 

shape of the distribution is obviously different. The sloping character 
of the top edge of the distribution is the result of neutron attenuation 
in the detector; the attenuation being of the order of e - * 84 . 
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Figure 18. Flight -time 
neutrons for a 6.01-cm. -thick c 
0 = 60°, and Y = 30°. Theresu' 
shown. 




I 




j 



PPARENT FLIGHT TIME, t; (nsec.) 

! 

distributions for 2-, 14.5-, and 100-MeV. 

letector. For these data L = 70 cm., 

1-ts for a. = 0.45 and 0.86 nsec, are 
J 

i 

I 

i 
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The results of these calculations are summarized further in Table IV. 
The detection efficiencies, obtained from the integral of the flight-time 
distribution over the range of flight times, are compared to the experi- 
mental values of Love et al ., 44 for a detector of approximately the same 
size. The differences in the values listed result primarily from the 
single scattering approximation used in this calculation to estimate the 
pulse -height spectrum. 

Improving the timing resolution for the configuration may be 
accomplished by increasing the flight path. For- example, for L = 200 cm., 
the overall timing resolution reduces by a factor of 2.7 which corresponds 
approximately to the increase in the flight time. There is, however, a 
corresponding decrease in the system efficiency arising from the reduction 
in solid angle subtended by the detector, and increased counting times 
would be required to gain comparable statistical accuracy in experimental 
data. 

It should also be noted that the values for the geometric resolu- 
tion given in Table IV differ. This arises since the jitter was sub- 
tracted in quadrature. This technique does not represent a precise 
way of eliminating this effect. 

II. COMPARISONS WITH CALCULATED FLIGHT-TIME UNCERTAINTIES 

The maximum relative flight-time uncertainties introduced by the 
target and detector dimensions were calculated in Chapter III under the 

44 T. A. Love et al.. Absolute Efficiency Measurements of NE-213 
Organic Phosphors for Detecting 1^.4- and 2. 6 -MeV. Neutrons, ORNL-3693; 
September 1966. 
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TABLE IV 

RESULTS FROM THE COMPUTATION OF 

c( VVV/[n T eff s p (v p ; V 0 ^ 

VERSUS t. 


(MeV. ) 

<V 

(nse' . ) 

* 

FWHM 
(nsec . ) 

GFWHM t 
(nsec . ) 

Total Time 
Resolution 
(*) 

Geometric 

Resolution 

(*) 



G j 

=0.43 nsec . 



2 

37.3 

2.51 

2.29 

6.7 

6.1 

14.5 

14.1 

1.19 

0.63 

8.5 

4.5 

100 

5.8 

1.02 

0.13 

17.5 

2.3 



( 7 . 

= 0.86 nsec. 





j 




2 

37-3 

2.92 

2.11 

7.8 

5.7 

14.5 

14.1 

2.11 

0.6l 

15.0 

4.3 

100 

5.8 

2.03 

0.15 

34.8 

2.5 



Computed 

Efficiency 

Zero Bias 


(MeV. ) 


Efficiency 

from Love* 

Efficiency** 


2 


28.5 

38.5 

56.9 


14.5 


18.1 

22.6 

27.8 


100 


6.5 


7-2 



* 

FWHM is full width at half maximum. 


t GFWHM is [(FWHMjP - (2.354 X 

J 

*T. A. Love et al., Absolute Efficiency Measurements of NE-213 
Organic Phosphors for Detecting 14,4- and 2.6-MeV. Neutrons, ORNL-3893, 
September 1966. 


•** 


These data are the same for a. = 0.43 and a. = 0.86. 

J J 
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assumption that all the dimensions have negligible effect except the one 
being evaluated. In order to approximate these same conditions in the 
computation of the physical flight-time distributions for a given dominant 
target or detector dimension, the remaining dimensions were made suffi- 
ciently small that their combined contribution to the flight -time 
distribution was small. The computer program logic prohibited their 
being set equal to zero. 

The numerical calculation of C(v ;v^, t^ ) is dependent on the 

existence of some fraction of instrumental time jitter. Recalling 

Equation IV-15, we note that o ., the standard deviation of the jitter 

J 

distribution, enters the calculation of the flight -time distribution 

principally through the evaluation of the two error functions. For ct. 

approaching zero, the difference between the error functions varies 

between 0 and 1 for combinations of quadrature points in the target and 

detector and the computed flight-time distribution "oscillates" with t. 

Since it is not possible then to use zero jitter, a . was made sufficiently 

J 

small that the full width of the jitter distribution was small compared 

to the full width at half maximum (FWHM) of the underlying "geometric" 

component of the flight-time distribution. The value of a ■ used in a 

given calculation was determined from estimates of the difference in the 

maximum and minimum flight times along a given direction f}. If the range 

of neutron flight times were such that the value of (t - t. )/a. and 

max -L j 

(t . - t.)/a. were greater than 4.215. the difference in the error 

min i J 

function was zero and the distribution goes to zero over this time 
region. Since the numerical calculation of the physical flight -time 
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distribution is sensitive to a , the optimum value was often determined 

J 

by observing the behavior of C(v ;v ,t.) as a function of a.. 

P ^ J 

Figure 1 9 shows the flight-time resolution, (At/t) induced by the 
detector geometry as a function of R^L and D/L. The flight -time resolu- 
tion was obtained from the ratio of the full width of half maximum (FWHM) 
to the calculated mean flight time, t , of the distribution; the instru- 
mental timing resolution component having been removed. 

Eliminating the effects of the jitter was accomplished in two ways. 

the width of the jitter distribution, given by 2.36 o was less than 

0 

10$ of the width of the combined jitter and geometric distributions, it 
was subtracted in quadrature. This method gives an imperfect result 
since the "geometric" component of the distribution is rectangular rather 
than Gaussian. 

The rectangular shaped distribution arises, in the case of the 
detector thickness, for the following reasons. In the absence of any 
neutron attenuation and with insignificant target dimensions, the 
probability for neutron detection given by Equation IV -h reduces to 

dt I(t,v n ,n) = Z(v n ) dj& D (t) . VI-1 

For ^(t) inconsistent with neutron detection in the detector, l(t,v n ,n) 
is zero. Otherwise, I is proportional to E(v n ) for all values of ^p(t) 
and the resulting distribution is flat across the range of allowed values 
of t. When the neutron attenuation cross section is introduced, the 
rectangular-shaped distributions reduces to a trapezoid as shown in 
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(O 2 10 - * '0° 

R 0 /L or d/l 


Figure 19 . The flight -time resolution induced by detector 
geometry versus R^/L and D/L for 2~, Ik *5-, and 100-MeV. 


81 


Figure 20. For this case, the detector thickness is 10 cm. and the 

attenuation is of the order of e -1 * 4 which produces the rapid drop off 

in the upper edge of the distribution. 

When 2.55 o. is small compared to the overall width of Cfv ;v ,t ) 
J p n 7 x' 

the error introduced by quadratic subtraction does not seriously affect 

the final answer; about 0.25$ for the case shown in Figure 20. For those 

cases where 10$ ^ 2.55 ct./FWHM ^ 57$^ the effects of the jitter were 

J 

removed using mathematical techniques proposed h y Kalisz 45 to determine 
the resolving time of a coincidence circuit from empirical coincidence 
curves. 

Provided the flight -time distribution is rectangular, we can write 

a i “ VI-2 


where 


= the standard deviation of the instrumental 
time jitter distribution, 

FWHM = the full width at half maximum of the computed 
flight-time distribution. 


Kalisz has tabulated the relationship between <j^ and t , where 2 (FWHM)t 
is the full width at half height of the true distribution corrected for 
the effects of smearing by a Gaussian distribution. 


45 Josef Kalisz, A Method for Determination of the Resolving Time 
and Efficiency of the Coincidence Circuit on the Basis of a Single 
Empirical Delayed Coincidence Curve, Institute of Nuclear Research, 
Warsaw, Report INR No. 652/IA/PL, September 1965 . 
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If the value of 2.35 a./FWHM > 35 1 °> the data was usually not 
considered valid and the case was re-evaluated using a smaller value for 

CT . . 

J 

We observe from Figure 19 that the computed values for the timing 
resolution, given by the data points, are in excellent agreement with 
the value predicted by the analytic expressions from Chapter III; the 
latter being given by the bold lines. For values of D/L > 5 x 10“ 2 we 
observe that the 2- and 14.5-MeV. points begin to drop away from the curve 
This apparent reduction in (At/t) is due chiefly to the method of computing 
the FWHM of the distribution. In all cases the FWHM was defined con- 
ventionally and, as can be seen from Figure 20, the width is small compared 
to the total width of the distribution. For the thick detector cases in 
which the shape is heavily influenced by the effects of neutron attenua- 
tion, observation of the actual predicted time distribution indicates 
no real discrepancy with the maximum flight-time distributions. 

The computed points for D/L fall slightly below the lines. This 
results since the lines were plotted using the equations derived in 
Chapter II in which the flight (At/t) is defined in terms of t = L/v 

n 

where L is the normal distance from the origin of the target coordinate 
system to the face of the detector. The computed values for (At/t) were 
obtained using the mean flight time. 

Also plotted in Figure 19 is the flight-time resolution as a 
function of R^/L. The agreement between the data is excellent in the 
region R^/L > 5 X 1CT 2 . Since the timing resolution is a function only 
of the radius and not the neutron energy, the points fall together as 
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expected. As in the case of the detector thickness, the flight-time 
distribution is almost rectangular. There is, however, a fall off in 
the top of the distribution which is inversely proportional to the 
flight time. This behavior in the shape of the distribution is pre- 
dicted by the analytic results when Equation IV-5 is integrated over 
all values of for a given flight time, t^. This integration is 
performed asuming the target dimensions, and detector thickness were 
negligible in extent. 

When Rjj/l < 6 x 10 -2 , the effects of the detector thickness used 
in the computation begin to dominate as shown by the open points. 

Figure 21 shows the flight -time resolution induced by the target 
thickness plotted as a function of T secY/L. The agreement between the 
computed values of (At/t), given by the points, and those predicted 
using the equation shown in the figure is favorable. (The reader is 
referred to Section II, Chapter III for the detailed calculation leading 
to this equation. ) As shown in the figure, some of the computed points 
have been corrected for the dominant detector contribution to the flight - 
time ie solution by subtracting these effects in quadrature. 

It is interesting to note the reversal in the expected flight - 
time resolution with increasing neutron energy. For this combination of 
target and detector angles (y = 30° and 9 = 60°), the extremum edges of 
the target can be used to approximate the maximum flight -time uncertainties 
for all neutron energies except for the particular case when cos 60° = 
v^/Vp. For this energy, the maximum flight time must be approximated 
using Equation III-15 etseq. 



FLIGHT-TIME RESOLUTION INDUCED BY TARGET THICKNESS (A///, %) 



'O 4 <0 3 10 2 10 * 10° 

T secifi/L 


Figure 21. 
versus T sec y/L. 
and Y = 30°. 


The flight-time resolution induced b y target thickness 
For 2.-, lk«5 and 100-MeV. neutrons with 0 = 60° 
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Figure 22 shows the shape of the computed distribution when 

T secY = 2.31 cm., L = 70-cm., and = 2 MeV. The time jitter comprises 

a significant portion of the total distribution. Since the distribution 

is rectangular, the effect of jitter may be removed using the Kalisz 

method. The FWHM is 0.433 nsec, and a. = 0.1 nsec. The resulting 

J 

"geometric" distribution obtained using the Kalisz method 0.432 nsec., 

while subtraction of < 7 . in quadrature yields 0.364 nsec.; about a 20 $ 

J 

variation. 

The comparison of the flight -time resolution is a function of the 
target radius requires a slightly different approach. In the derivation 
of the maximum flight -time uncertainties arising from the target radius, 
we assumed uniform beam density over the target radius. However, in 
deriving C(v ;v n ,t^ ), we assumed that, in the absence of multiple scat- 
tering of the protons, the target radius corresponded to the radial 
dimension of the beam. The beam density was approximated by a Gaussian 
distribution in the x-z plane. In the numerical integration leading to 
C(v ;v n ,t^), the limits for the y integral were defined by y = ± 2 . 5 ( 7 , 
where <7 is the standard deviation of the beam distribution. 

To demonstrate the effects on the flight -time distribution of the 
target radius, we must introduce a different, but consistent, notation 
into the equations for leading to (At/t) given in Chapter III. 

•I'm 

The "radius" of the target was approximated by the FWHM of the 
beam distribution. Considering the rotation angle Y, the effective 
target radius was defined by 

( R T ) e ff = 1,18ct S6CY = °y ' 


VI-3 
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49 50 51 52 

JITTER SMEARED APPARENT FLIGHT TIME, t, (nsec) 

Figure 22. The physical flight-time distribution for the case 
where T(|) = 2.0 cm. and all other target and detector parameters have 
been minimized. For these data F = 2.0 MeV., L = 100 cm., 0 = 60°, 

Y = 30 , and a. = 0.1 nsec. 
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Equation III-20 may now be written with a' replacing R.^. 

Figure 23 shows the flight -time resolution as a function of a* / L. 
The points give the computed values for (At/t) obtained from the flight - 
time distribution while the lines were obtained using Equation III-20. 
For these cases, the instrumental jitter was subtracted in quadrature 
since the underlying "geometric" distribution is Gaussian. 

In the numerical computation leading to Figure 23 the integration 
over the p coordinate of the detector was performed on the rear face 
of the detector. For the small detector used in these studies, the 
uncertainty in the final results due to the edge effects is negligible. 

III. THE EFFECTS OF TIME SLEWING AND 
THE TIME-CHANNEL RESPONSE FUNCTION 

Figure 24 shows the effects of the time slewing on the physical 
flight-time distribution for 14.5-MeV. neutrons obtained in Section I 
of this chapter. These data were obtained from Equation IV-38 using 
the time-slewing response function for 14.5-MeV. neutrons given in 
Figure 15 . 

Comparing the flight -time distribution from Figure l8 with these 
data, we observe significant broadening the original distribution in 
time with a corresponding truncation of the amplitude. The integral 
over the measured flight time t yields the efficiency of the detector 


obtained in Table IV. 



FLIGHT-TIME RESOLUTION INDUCED BY THE "EFFECTIVE" TARGET RADIUS {A t/t } %) 
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10 1 2 5 10 2 2 5 ICf* 

o^/L 


Figure 23. The flight-time resolution induced hy the "effective" 
target radius versus cr'/L f° r 2-, 14.5-, and l(X)-MeV. neutrons. 

9 = 60° and Y = 30°. y 



15 


MEASURED FLIGHT TIME, t m (nsec.) 

Figure 2k. The probability per nsec, per incident proton per 
MeV. of neutron energy that a neutron has a measured flight time at t 
versus t . the measured flight time. 
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Figure 25 shows the time-channel response functions obtained from 
the slewing -smeared distribution in Figure 2k. These response functions 
are obtained from Equation IV-37 using an analyzer calibration of 
0.500 nsec ./channel. The values obtained in this analysis are compared 
with the experimental points of Love et al. 4 6 for approx ima tely 14-MeV. 
neutrons at a flight path of 40 cm. The experimental data were 
normalized as well as possible for the differences in the flight path 
and the analyzer .calibration. 


46 


Love, op. cit., p. 22. 
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Figure 25. The time-channel response functions for 14.5-MeV. 
neutrons . 
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APPENDICES 


APPENDIX A 


CALCULATION OF t : THE INTERVAL FROM TIME ZERO TO THE 

P 

(p,n) INTERACTION IN THE TARGET 


The time spent by the proton in the target is given by the relation 


t 

P 




+ t 

e 


A-l 


where 

dx = the incremental distance travelled by the proton in 
the target, 

v (x) = the proton velocity at the distance x in the target, 

P and 

t = the time required for the proton to travel from the 
plane passing through x = 0 to the entrance face of 
the target at velocity v . (t is usually a negative 
number.) ^ 

The proton loses energy in the target through numerous atomic 
collisions. Its flight time through the target is a function of the 
incremental energy loss per cm., dE/dx. In the nonrelativistic approx- 
imation, the proton energy is given by the familiar equation E = 

(l/2) mv^. Differentiating with respect to x we obtain 

(-dE/dx) = mVp (-dv^/dx) . A-2 
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(-dE/dx) is a positive number; the minus implying a decrease in energy 
with increasing x. Expanding v“ 1 (x) in a Taylor series gives 

Vp 1 ^) = v p 1 (°) " ( v p 2 x ) (dv p /dx) + ... . A-3 

Substituting dv /dx from Equation A-2 into A-3 and introducing 
the result into Equation A-l gives 

\ = J* P F" 1- ~ ( -dE/dx ) J p xdx + t . A-k 

* 0 p mv 3 0 

P 

JL is the proton penetration depth prior to the (p-n) interaction given 
"by Equation II-3. Noting that mv^ = and completing the integration, 

we obtain the result 


l i, 2 

\ ^ + WW + t • A-5 

p p p 

The convention was adopted that the zero of time for a proton is 
at the instant it would have passed through the plane x = 0 if the target 
were absent. Then, 


t 


e 


x 



v 

P 


A -6 


100 


Substituting this result into Equation A-5 yields 


‘p 25 f + ( - dE/ax) 

P P P 


A-7 



APPENDIX B 


CALCULATION OP (1' • h), A , AND JL^ 

The quantities (i' • Q), and £ T introduced in Chapter IV in the 
analysis leading to the physical flight-time distribution of detected 
neutrons will be derived here in terms of the known quantities L, x, y, 
y ) z , Q } 1, and t 2 • 

A 

I. THE DERIVATION OF A and (i' • fi) 

The distance travelled by the neutron from its birth point in 
the target to the face of the detector is JL Introducing Sc and y? 
the target interaction coordinates in a coordinate system rotated 
through an angle © in the x-y plane, 

x = x cos© + y sin© 

* "R 1 

y = y cos 0 - x sin© 

we can write the cartesian components of Q in the rotated coordinate 
system as 


Q / = L - x 
x 

= y' - y ; b-2 

n , = z' 

z 
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if we make the identification 


o — • 


It follows from the Pythagorean Theorem that 


l t = J (L - X? + (y' - y f * ( Z ' f . 


f f \ A 

(i • n ; is the cosine of the direction angle of fi relative to the x 
coordinate system. From B-2 , 



It can be shown from similar arguments 
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II. DERIVATION OF the PATHLENGTH OF THE 
SECONDARY NEUTRON IN THE TARGET 

The length of & is dependent on the target -detector orientation 
and separate calculations are required for the target in transmission, 

|© - y| ^ rt/2, and in reflection |© - y| > jt/2. 

The Target in Transmission 

Neutrons born in the target at P T are emitted at an angle 0 with 
respect to the proton beam axis. JL T is obtained from the normal equation 
of the plane of the exit face of the target. According to this equation, 

Ji T coscp = |on| B-6 

where ON is the normal from P T to the exit face and cp is the angle 

^ A 

between ON and Cl . The length of ON is 

|on| = t 2 - x COSY - y sinY • B-7 

coscp is obtained from the equation 

coscp = COSQ^ COSQ^ + COS^ COSf^ + COS COS 7 2 £_8 


where (c^ ^ y^ ) are the direction angles of and (a, ^ y ■ ) those of ON; 
all angles being measured relative to (x'y'z ' ) along the direction 6 


its direction cosines are given by Equations B-4 and B-5. The direction 
cosines for ON are given by 


cosa 

i 

/\ 
= ON 

A 
. / 

. 1 

= cos(Y 

- 0) 

cosP i 

\o 

ii 

* J 

= sin(Y 

- 0) 

cosy^ 

= oil 

• k" 

= 0 



where ON" is the unit vector along ON. 'ON * k' = 0 since ON is constrained 
to the x-y plane. Then Equation B -6 can be written 


(t 2 secY -x-y tanY ) cosY 




B-10 


/• ^ 


r ~ x cos(e - y) - - y sin(0 -y) 

Mr S 


which gives the magnitude of Si ^ for the target in transmission. 

Target in Reflection 

When 1 0 - Y | > it/ 2, the derivation for & ^ follows the same 
procedure . In this case, |on| is given by cosY where & is given by 
Equation II-3 . Then, 


•^ip 


(x + y tanY - t x secY ) cosY 
cos ( 0 - Y) - y sin(0 - Y) 


B-ll 


for the target in reflection 
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